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Abstract 

A  constitutive  equation  is  proposed  with  a  view  to  describing  the  rate 
dependent  mechanical  response  of  metals  at  high  temperatures.  The  equation 
is  of  the  endochronic  type  and  derives  its  physical  foundations  from 
deformation  kinetics.  Of  importance  is  the  fact  that  hardening  is 
associated  with  a  change  in  the  energy  barriers  brought  about  by  the 
inelastic  deformation  of  a  metal.  The  equation  is  used  to  describe  the 
results,  by  Ohno  and  his  associates,  of  experiments  on  the  creep  response 
of  metals  to  piece-wise  constant  stress  histories.  The  metal  in  the 

yO- 

present  case  is  304  stainless  steel  at  600°^ C.  It  is  shown  that  the  theory 
gives  analytical  results  that  are  in  close  agreement  with  experiment. 

Of  consequence  is  the  fact  that  the  constitutive  equation  applies  to 
three-dimensional  stress  or  strain  histories  and  is  thus  not  limited  to 
those  stress  histories  associated  with  creep. 


1.  Introduction 


In  the  present  paper  we  develop  an  endochronic  theory  of  viscoplasticity 
which  accounts  for  the  history  of  strain  and  strain  rate  on  the  stress  response 
of  metals  at  high  temperatures.  The  theory  is  based  on  the  concepts  of 
endochronic  plasticity  (see  typically  Ref.'s  [1],  [2],  and  [3]),  however,  the 
increment  of  intrinsic  time  scale  is  no  longer  proportional  to  the  plastic 
strain  path  but  depends  also  on  the  rate  at  which  the  path  is  traversed.  The 
development  of  the  theory  is  dealt  with  at  length  in  the  subsequent  sections. 

The  resulting  constitutive  equation  is  used  to  analyze  the  creep  response 
of  304  stainless  steel  to  piece-wise  constant  shear  stress  histroies  at  600°C 
and  to  compare  the  results  with  the  experimentally  determined  creep  response  of 
the  same  material  at  this  temperature  as  reported  by  Ohno  et  als.  in  Ref  [4], 

The  constitutive  behavior  of  metals  at  high  temperature  where  the  strain 
rate  sensitivity  of  the  mechanical  response  cannot  be  ignored,  has  been  the  sub¬ 
ject  of  extensive  research  in  recent  years.  We  do  not  wish  to  give  in  this 
paper  an  exhaustive  review  of  the  literature  on  this  subject  but  merely  cite 
references  which  are  typical  of  the  enormous  amount  of  work  which  is  being  done 
in  this  field.  In  this  context,  the  works  of  Chaboche  [5],  Krieg  [6],  Malvern 
[7],  Haisler  [8],  Bradley  [9],  Leckie  [10],  Krempl  [11],  Walker  [12],  Miller 
[13],  and  Hart  [14],  among  others  are  mentioned.  We  also  wish  to  cite  the  works 
of  Ohashi  et  als.  [15],  and  Murakami  and  Ohno  [16],  and  Ohno  et  als.  [4]  who 
have  been  enjoying  a  measure  of  success  in  using  a  creep  hardening  surface 
theory  in  describing  the  creep  response  of  metals  to  piece-wise  constant  stress 
histories,  a  subject  with  which  we  will  be  dealing  in  this  paper.  This  aspect 
of  the  mechanical  response  of  metals  has  given  rise  to  greater  difficulties 
than  say,  the  stress  response  to  piece-wise  constant  strain  rate  histories. 
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The  work  of  Krausz  [17]  and  his  co-workers  also  occupies  a  significant  position 
in  the  literature  because  of  its  more  fundamental  nature  in  that  a  microscopic 
theory  of  deformation  kinetics  is  used  to  gain  understanding  of  the  mechanical 
response  of  metals  at  the  bulk  level. 

In  our  initial  approach  to  the  subject,  with  specific  attention  to  the 
strain  (creep)  response  to  piece-wise  constant,  pure  shear  stress  histories,  we 
used  a  strictly  phenomenological  approach,  in  the  context  of  an  endochronic 
theory.  Specifically  in  one  dimension  we  probed  the  data  with  the  constitutive 
equation 


ep  =  /  J(z-z' )  ^|t-  dz' 


(1.1) 


Henceforth  convolution  integrals  such  as  the  one  appearing  on  the  right  hand 
side  of  eq.  (1.1)  will  be  given  in  symbolic  form  according  to  eq.  (1.1a): 


z  a.  def 

/  J(z  -  z' )  dz'  =  J(z)*ds 
0- 


(1.1a) 


D  C 

In  equation  (1.1)  eK  has  the  same  connotation  as  e  ,  i.e.,  it  is  the  inelastic 


shear  strain  where 


•p  *c  s 

e  =  e  =  e  -  ■=— 

2Ur 


(1.2) 


and  uQ  is  the  elastic  shear  modulus.  We  caution,  however,  that  other  defini¬ 
tions  of  ec  have  been  used  in  the  literature.  The  intrinsic  time  z  was  defined 
by  equation  (1.3)  where 


f(c,c) 


(1.3) 


:c,  , 


and  c  =  [e  |  is  the  usual  fashion.  The  dependence  of  f  on  c  lends  "strain  rate 


sensitivity"  to  the  equation,  which  otherwise  would  be  strain  rate  insensitive. 


While  equations  of  the  type  (1.1),  (1.2)  and  (1.3)  were  shown  to  give  satis¬ 
factory  results  in  the  case  of  constant  strian  rate  histories  as  demonstrated  by 
Wu  and  Yip  [18,19]  and  Lin  and  Wu  [20,21],  we  found  that  these  equations  did  not 
prove  satisfactory  in  the  case  of  piece-wise  constant  stress  histories.  With 
specific  reference  to  the  data  of  Ref.  [4],  it  was  found  that  f  had  to  be  essen¬ 
tially  independent  of  z,  to  account  for  the  periodic  creep  response  to  piece- 
wise  constant  cyclic  histories.  Thus,  limiting  f  to  a  dependence  on  c  only 
resulted  in  a  gross  overestimate  of  the  creep  strain  under  cyclic  conditions. 
Correcting  f  so  as  to  match  the  data  gave  rise  to  oscillations  in  f  which  could 
not  be  accounted  for  by  means  of  equation  (1.3).  Furthermore,  the 
phenomenological  approach  did  not  give  any  hint  as  to  the  physical  mechanism(s) 
responsible  for  such  fluctuations  in  f.  To  overcome  the  difficulties  we  appealed 
to  the  theory  of  deformation  kinetics  in  the  context  of  the  internal  variable 
theory.  The  latter  is  treated  briefly  in  Section  2  while  the  former  is  repre¬ 
sented  in  detail  in  Section  3. 


2.  Internal  Variable  Theory 


Irreversible  thermodynamics  of  internal  variables  is  now  a  well  established 
field  so  we  proceed  to  give  a  very  breif  outline  of  the  theory  for  the  sake  of 
completeness.  We  limit  ourselves  to  small  strain  fields.  In  the  Helmholtz  for¬ 
mulation  the  thermodynamic  state  is  described  by  the  free  energy  density  <f>  which 
is  a  function  of  the  strain  tensor  £  the  temperature  T  and  n  interval  variables 

r 

q  which,  for  thermomechanical  processes,  are  tensors  of  the  second  order. 

w 

The  stress  <3  and  the  entrpy  density  ^  are  then  given  by  equations  (2.1)  and 

(2.2) 


1  = 


9\J» 

FT 


(2.2) 


In  the  case  of  a  spatially  uniform  thermal  field  the  condition  of  positive  rate  of 
irreversible  entropy  gives  rise  to  the  inequality 


3il>  •  r  ^  •  r 

-  •  q  >  0  ,  iq  l  *  0  (2.3) 

9qr 


r,  not  summed. 

In  the  linear  version  of  the  endochronic  theory  4>  is  a  quadratic  function  of 

r 

its  arguments  and  the  evolution  equations  for  q  have  the  form 


3ij> 

9q 


-  +  b 
r  „r 


dqr 

~dz 


(2.4) 


r  not  summed,  where  n  such  equations  exist,  one  for  each  variable  q  ,  and  b 

~  ts  ^ 

are  positive  definite  tensors  of  the  fourth  order. 


r  ✓ 


,»  ,*  V  V*  /  V*  v  v  v 


i 
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When  isothermal  conditions  prevail  and  the  solid  is  initially  isotropic  and 


the  tensors  ty  are  constant,  then  equations  (2.1)  and  (2.4)  combive  to  give  rise 


the  integral  constitutive  equations 


s  =  2u(z)*de 


(2.5) 


0  =  3K(z)*de  (2.6) 

where  s  is  the  stress  deviator,  e  the  strain  deviator  and  0  and  e  the  hydro- 
estatic  stress  and  strain  respectively.  The  kernels  p(z)  and  K(z)  are  sums  of 
positive  decaying  exponential  functions,  i.e.. 


-arz 


fMz)  =  l  e  ,  K(z)  =  l  K  e 


-Xrz 


(2.7a ,b) 


r  r 

where  U  ,  a  ,  K  and  X  are  all  non-negative, 
j  r  r  r  r 

In  the  generalized  endochronic  internal  variable  theory  the  evolution 
equations  are  expressed  in  terms  of  the  intrinsic  times  of  the  mechanisms  of 


internal  motion.  See  Ref  [22].  Specifically  to  each  qr  the  theory  ascribes  an 

tf 
dqr 


intrinsic  time  zp  such  that  the  equations  of  evolution  become 


+  b  • 


•  — —  =  0  (r  not  summed) 


(2.8) 


r  =  1 ,  2  . . .  n. 

In  the  endochronic  theory  of  plasticity  of  metals  as  it  has  been  used  in  the 

past 


Z1  "  Z2  "  •”  =  Zn  =  Z 


(2.9) 


and  the  elastic  bulk  modulus  K  in  constant,  i.e,  the  material  is  plastically 
incompressible.  Also, 


dz  =  ,  dc  =  ude13! 


fk) 


(2. 10a ,b) 
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* 


r 

[ 


I 

t 


and 

P  d-§ 

dep  =  de  -  p— -  (2.11) 

y0 

where  u0  is  the  elastic  shear  modulus.  The  function  f  is  positive  and  non 
decreasing.  Thus  in  the  case  of  plasticity,  and  in  the  context  of  the  above 
assumptions,  equations  (2.5)  and  (2.6)  have  the  form 


s  =  2v>0(z)*de  (2.12) 

G  =  3Ke  (2.13) 

Substitution  of  equation  (2.11)  in  equation  (2.12)  gives  an  equivalent  consti¬ 
tutive  equation  which  relates  s  directly  to  the  history  of  eP.  Thus 

s  =  2p(z)*deP  (2.14) 


See  Valanis,  Ref.  [23].  The  relation  between  p  and  jju  is  given  in  terms  of  their 
Laplace  transforms  in  equation  (2.15). 

p"  (1  -  —  )  =  IT  (2.15) 

p0 

It  has  been  found  that  in  the  case  of  metals  at  room  temperature  p  and  f  are 
well  represented  by  the  relations 


1  ~kz  r  -i  -bz  . .  ,,  .  . 

p  =  p0  ~  e  .  f  =  1  -  ae  (2.16a,b) 

where  p ^ ,  a,  a,  b  are  positive  and  k  is  non-negative. 

In  the  Gibbs  formulation  the  thermodynamic  state  is  described  by  the  free 
energy  density  <fr,  which  is  a  function  of  the  stress  tensor  a  the  temperature  T 

p 

and  n  internal  variables  q  which,  again,  are  second  order  tensors.  The  func¬ 
tion  (j>  is  related  to  by  the  equation 


(2.17) 


The  counterparts  of  equations  (2.1)  and  (2.2)  are 
1  =  '  ^ 


and 


3<j> 

Tf 


(2.18) 


whereas  the  positive  rate  of  irreversible  entropy  gives  rise  to  the  inequality 


3<t>  *  r  ■  r 

— —  •  q  >  0 ,  !!  q  il  *  0 

3qr  ~ 


(2.19) 


3  d) 

for  each  r,  where  -  - —  is  the  internal  microforce  for  the  internal  mechanism  r 

3qr 


while  the  evolution  equations  for  q  are,  in  the  case  of  the  generalized 
endochronic  formulation 

r 


3  <(> 


r  +  br 

a3  “  r 


dor 

dz 


0 


(2.20) 


If  we  now  stipulate  that  4>  is  quadratic  in  its  variables  and  dzr  =  dz  for 
all  r,  and  b  are  constant  then  in  the  case  of  isotropic  materials  and  isothermal 
conditions,  equations  (2.17)  and  (2.20)  combine  to  give  the  constitutive 
equations 

e  =  \  L(z)*ds  (2.21 ) 

g  =  j  M(z)*do  (2.22) 

When  plastic  incompressibility  applies  N  is  a  constant  and  equal  to  (see 
equation  (2.13).  Also  u  and  L  are  related  by  equation  (2.23). 


8 


u (z)*dL  =  H(z) 

where  H(z)  is  the  unit  step  function. 
In  view  of  equation  (2.11) 


p  1  ^ 

e  7 T~ 

A#  ri  £  Hg 


(2.23) 


(2.24) 
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Thus 


r*  =  7  >Hz)*dS 


(2.25) 


where 


J(z) 


2v, 


H(z)  +  L(z) 


(2.26) 


[i 

I 


,V 

,% 

A 


The  functions  J(z)  and  p(z)  are  also  related  as  shown  in  equation  (2.27) 


J (z)*dp  =  H ( z )  (2.27) 


Spectrum  of  Intrinsic  Times 


Of  the  n  internal  mechanisms  let  group  r  have  an  intrinsic  time  zr>  Then  to 

the  group  r  of  internal  variables  n  in  number,  there  will  correspond  an  intrin- 

m  r 

sic  time  z  .  Evidently  £  n  =  n  where  m  is  the  number  of  groups. 
r  r=l  r 

A  straightforward  analysis  using  equations  (2.16)  (2.17)  and  (2.20)  when  <f>  is 

( r ) 

a  quadratic  isotropic  function  of  e..  and  q..  and  b  are  constant  isotropic 

TJ  iJ  .r 

tensors  leads  to  the  equations 
„  n 


eH  =  l  J  (z  )*ds 
~  r=l  r  r  " 


(2.28) 


a 


I 


l\| 

is] 


1 


where  plastic  incompressibility  has  been  observed.  They  physical  interpretation 


of  equation  (2.28)  is  that  each  of  the  m  groups  of  mechanisms,  say  r,  contributes 


to  the  total  strain  a  partial  strain  ej.  such  that 


I  eP=ep 

r=l  ~  ~ 


(2.30) 


where 


er  =  Jr(zr)*ds 


(2.31) 


The  intrinsic  time  zr  is  related  to  e  by  the  equation 


where  f  is  the  hardening  function  of  graph  r. 


The  need  for  this  more  general  approach  which  was  presented  in  a  previous 


reference  [22],  has  been  discussed  in  the  introduction  and  ha§  to  do  with  the 


fact  that  one  intrinsic  time  is  just  not  sufficient  to  describe  the  creep 


response  of  metals  to  piece-wise  constant  stress  histories.  The  physical  justi¬ 


fication  for  this  possibility  is  discussed  in  the  section  on  deformation  kine¬ 


tics,  but  simply  put,  it  means  that  the  hardening  function  f  is  not  the  same  for 


all  mechanisms  at  high  temperatures,  though  the  assumption  of  an  f  common  to  all 


q  suffices  at  room  temperatures  as  demonstrated  in  our  work  on  endochronic 


plasticity.  The  appeal  to  deformation  kinetics  is  necessitated  by  the  desire  to 


determine  how  f  is  influenced  by  the  micro-mechanical  process  which  accompanies 


the  inelastic  deformation. 


The  theory  of  deformation  kinetics  was  founded  by  Eyring  [24]  and  its  appli¬ 
cation  to  rate  processes  has  been  pursued  by  Eyring,  Krausz  [24]  and  their 
co-workers  with  a  great  deal  of  success. 

In  keeping  with  the  ideas  of  deformation  kinetics  we  attribute  macromotion 
to  additive  effects  of  micromotions  brought  about  by  local  distortion  of  atomic 
"energy  barriers".  At  this  point  it  is  essential  that  we  distinguish  between 
diffusion  of  particles  and  diffusion  of  dislocations  or  vacancies.  The  distinc¬ 
tion  may  be  stated  most  simply  in  terms  of  the  mean  free  path  £  of  a  particle. 

In  the  case  of  particle  diffusion  £  is  large  compared  to  its  counterpart  in  the 
case  of  dislocation  or  flaw  diffusion.  If  a  particle  travels  n  units  of  distance 
"a"  before  coming  to  rest,  then  £  =  na.  However,  if  a  flaw  travels  n  units  of 
distance  a,  the  mean  free  path  of  a  particle  is  still  a,  because  a  different 
particle  partakes  each  time  in  the  motion  of  the  flaw.  Thus  in  the  case  of  par¬ 
ticle  diffusion  n  >>  a.  Most  important,  however,  is  the  fact  that,  in  either 
case,  each  unit  of  motion  consists  of  an  atom  moving  across  an  energy  barrier. 

In  addition  we  would  expect  that  in  the  case  of  flaw  or  dislocation  diffusion 
the  energy  barriers  would  be  lower  than  those  of  particle  diffusion.  In  fact, 
the  activation  energy  of  self  diffusion  is  lower  at  low  homologous  temperatures 
(where  dislocation  motion  is  dominant)  than  at  higher  temperatures  where  par¬ 
ticle  diffusion  dominates  the  process. 

To  apply  the  ideas  of  deformation  kinetics  to  the  viscoplastic  deformation 
and  flow  of  solids  we  appeal  to  a  simple  atomic  model  whereby  prior  to  the 
application  of  stress  each  atom  of  the  solid  is  situated  at  the  bottom  of  a  sym¬ 
metric  potential  well.  A  typical  potential  well  with  the  accompanying  local 
potential  surface  is  shown  in  Figure  1  by  a  solid  line.  The  forward  and  back- 

p 

ward  barriers  are  equal  and  both  have  a  height  . 


When  stress  is  applied, 


the  atoms  will  be  displaced  from  their  initial  positions  and  the  local  potential 

surface  of  an  atom  will  distort.  The  distorted  potential  surface  is  also  shown. 

.  r 

The  effect  of  the  distortion  is  to  reduce  the  forward  barrier  by  an  amount  w^ 

,  r 

and  increase  the  backward  barrier  by  an  amount  w^  .  It  will  be  shown  that  this 
type  of  barrier  distortion  will  give  rise  to  an  average  forward  motion  of  the 

,  r 

atoms  occupying  potential  wells  with  barriers  £q. 

Boltzmann  Statistics.  To  determine  quantitatively  the  effect  of  barrier 
distortion  on  the  mean  atomic  motion  we  appeal  to  Boltzmann  statistics. 
Accordingly  the  probability  of  finding  an  atom  in  an  energy  state  e.  is  given  by 
equation  (3.1): 

-8e.j 

Pi  =  «e  (3.1) 

where 

"Bei  1 

a  =  1/1  e  ,  8  =  yt  (3.2a,b) 

i 

in  the  usual  notation. 

Let  Nr  be  the  number  of  atoms  occupying  potential  wells  with  initial  energy 

r 

barriers  eQ  .  The  probability  that  an  atom  is  in  an  energy  state  greater  than 

r  .  r  . 
is  Pq  where 

pJJ  =  l  <*  exp  (-BeJ)  (3.3) 

r  r 

ei>£0 

States  e!"  such  that  eT  >  we  have  called  activated  states  [25]  differing  from 

Eyring.  Thus,  the  number  of  atoms  in  an  activated  state  is  N  p!T  .  As  the 

r  u 

barriers  are  symmetric  the  probability  that  an  atom  will  move  forward  is  equal 

to  the  probability  that  it  will  move  backwards  so  that  the  net  motion  (average 

displacement)  of  the  atoms  N  is  zero. 

r 


vl-'/.  vlvi 


When  the  potential  energy  surface  is  distorted  the  number  of  atoms  Ar  par¬ 
taking  in  a  forward  motion  is  now  changed  to 


A  =  N  I  a  exp  (-eeT) 

r  r  r 

ei>e0_Wf 

while  the  number  Br  of  atoms  partaking  in  as  backward  motion  is 


(3.4) 


B  =  N  l  a  exp  (-8eV) 

r  r  r 

e  .  >e  +w . 

1  0  b 


(3.5) 


The  net  number  of  atoms  that  partake  in  a  forward  motion  is,  thus,  A  -  B^  where 


r  r  r 
e  .  <e  +wL 
i  0  b 

A  -  B  =  N  l  ct  exp  (-8eJ) 
r  rr  r  r 

e ->ek"w* 
i  b  f 


(3.6) 


To  evaluate  the  sum  on  the  right  had  side  of  equation  (3.6)  we  shall  assume  that 
wj  and  wbr  are  both  small.  In  this  case  we  represent  the  distribution  of 

energies  eT  in  the  vicinity  of  in  terms  of  the  local  tangent  to  the  distribu- 

•  r  .  . 

tion  at  eQ  by  writing 


Ei  =  eo  +  kl"(i  -  i0) 


(3.7) 


where  i'q  is  the  value  of  i  at  e.  =  and  kr  is  the  slope  of  the  distribute 
which  is  a  function  (in  general)  of  .  See  Figure  2. 

Substitution  of  the  relation  (3.7)  in  the  sum  on  the  right  hand  side  of 
equation  (3.6)  leads  to  the  simple  expression 

o.r 


-Be’ 
e  0 


A  -  B  =  2N  a 
r  r  r 


sin3„w  _ 
r  r 


(3.8) 


1  -  e 


where 


(3.9,10) 


vJ 


r  wfwb 


r  wb+  wf 

Aeo 


r  r 

£  “3k 

We  note  that  the  terms  e  0  and  1  -  e  are  related  to  the  initial  state,  i.e., 
the  barrier  height  and  the  energy  distribution,  while  the  terms  sin  Bwr  and 

p 

exp  (-0A£q)  are  brought  about  by  barrier  distortion. 

The  mean  velocity  vP  relative  to  the  lattice  of  the  atoms  in  group  r  may  now 

p 

be  calculated  in  terms  of  the  barrier  distortion  parameters  and  A eQ  .  If 
X  is  the  mean  lattice  distance  and  tr  is  the  average  time  taken  by  the  atoms  of 

p 

group  r  to  traverse  that  distance  across  the  barrier  eQ  then 


vp  =  (X/r  HA  -  B  )/N 
r  r  r  r  r 


(3.11) 


We  now  define  an  internal  variable  q^  by  the  relation 


%  =  vp/x 

r  r 


(3.12) 


Evidently  in  view  of  equation  (3.11) 


q  =  (A  -  B  )/N  t 
r  r  r  r  r 

In  view  of  equations  (3.8)  and  (3.13) 
.r 

“PS 

2a  e 

q 


(3.13) 


-S£0  -8AeS 

- =-  (— )e  si nh  6  w 

-Bkr  Tr  r 


(3.14) 


1  -  e 

In  so  far  as  steady  creep  is  concerned  the  assumption  is  usually  made  that 

of  all  the  operating  mechanisms  only  one  survives  in  the  steady  state,  i.e., 

r  =  1  and  w^  is  proportional  to  the  stress  (in  one-dimensional  stress  fields). 

However,  in  the  case  of  transient  creep  it  is  the  local  microforce,  i.e., 

3<j> 

-  on  the  group  r  that  will  determine  the  barrier  distortion.  Thus,  following 
Ref.  [25],  we  let  wr  be  proportional  to  -  according  to  equation  (3.15) 


motion  of  the  group  r  of  paricles  facing  a  potential  barrier  g  . 

Discussion  of  equation  (3.17).  As  we  pointed  out  equation  (3.17) 
establishes  a  physical  meaning  for  the  internal  variables  in  that  qr  is  the  mean 
displacement  relative  to  the  lattice  of  a  group  r  of  particles  facing  a  poten- 

r 

tial  barrier  of  magnitude  eQ  .  The  above  equation  was  published  by  Valanis  and 

r  r 

Lalwani  in  Ref.  [25],  with  AeQ  =  0.  The  appearance  of  the  term  AeQ  in  equation 
(3.1)  was  inferred  as  a  result  of  our  effort  to  describe  analytically  the  creep 
response  of  306  stainless  steel  to  piece-wise  constant  stress  histories.  This 
will  be  discussed  in  Section  4. 


The  time  to  traverse  the  barrier,  i.e.,  is  also  of  central  importance  in 
equation  (3.17).  Eyring  used  simplifying  assumptions  to  arrive  at  the  conclusion 
that  x  is  proportional  to  the  square  root  of  the  ambient  temperature.  However, 
one  can  ,how  that  it  depends  at  least  in  part  on  the  barrier  shape  and  height 
(Ref.  [25]).  In  this  work  we  have  found  that  is  also  sensitive  to  the  plastic 
strain  rate.  This  is  to  be  expected  sind  Tr  depends  on  the  barrier  confor¬ 
mation,  which  in  turn  depends  on  the  plastic  strain.  The  rate  of  plastic  strain 
affects  the  rate  of  barrier  distortion  which  must  affect  the  traversal  time  t  . 


15 


Consider  now  two  processes  (a)  and  (b)  the  first  of  which  is  proceeding  at  a 
faster  plastic  strain  rate  than  the  second.  With  regard  to  the  forward  motion 
of  a  particle,  the  height  of  the  barrier  will  be  diminishing  faster  in  case  (a) 
than  in  (b),  so  that  the  forward  moving  particle  will  be  encountering  a  con¬ 
sistently  lower  barrier  in  case  (a)  than  in  (b).  It  follows  that  the  time  to 
cross  the  barrier  in  case  (a)  will  be  shorter  than  in  case  (b).  Thus 


t  <  tu  whenever  c  >  5. 
a  b  a  b 

The  above  inequality  will  be  satisfied  if 


g(5) 


(3.18) 


where  tq  is  a  constant  and  g  is  a  monotonically  increasing  function  of  c.  In 
this  work  we  have  set 


9(?)  =  r,.m  (3.18a) 

c 

where  m  is  a  material  constant. 

Deformation  kinetics  is  brought  into  accord  with  linear  irreversible  ther¬ 
modynamics  if  in  equation  (3.17)  the  argument  of  the  hyperbolic  sine  is  suf¬ 
ficiently  small  for  the  approximation 


,  36  ,  36 

si nh  (8  C  T—  )  ~  6  C 


>  3q„  '  B  r  3q. 


(3.19) 


r  r 

to  be  appropriate.  In  this  event  equation  (3.17)  becomes  a  standard  linear  evo¬ 
lution  equation,  i .e. , 


br,r  ♦  |A-  =  0 


(3.20) 


where 


b  r  JL.  e 

u  r\  r»  nr* 


8  Ae, 


00  8C 


(3.21) 
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In  view  of  equations  (3.18a),  (3.20)  and  (3.21)  one  finds  that  the  "endochronic" 
form  of  equation  (3.20)  is 


(3.22) 


where 


r  .  r  T0 

!o  "  oo  sc 

r 


(3.23) 


and 


.  _  dz 

Zr  "  cmf 


(3.24) 


f  =  e 
r 


(3.25) 


Thus  in  deformation  kinetics  terms  the  rate  sensitivity  is  attributable  to  the 
time  to  cross  the  barrier  while  the  hardening  (softening)  is  related  to  the 


change  in  the  mean  height  of  the  barrier  as  a  result  of  the  stress  history. 

p 

Thus  if  Acq  increases  the  material  hardens  while  if  it  decreases  the  material 
softens  in  accord  with  our  physical  intuition  regarding  such  processes. 


4.  Analysis  of  Piece-wise  Constant  Stress  Histories  in  Pure  Shear. 


Wo  begin  with  the  integral 


eP  =  J(z)*ds  (4.1 ) 

where  eP  represents  a  shear  creep  strain  component,  s  is  the  corresponding  shear 
stress  component  and  J  the  appropriate  creep  function.  As  usual 


dz 


d£ 

g(S)f 


(4.2) 


where  g  is  the  rate  sensitivity  function  and  f  the  hardening  function.  Also 


dc  =  k|deP|  (4.3) 

where  k  is  an  appropriate  scalar  constant.  Typically,  it  eP  denotes  a  creep 
shear  strain  component  and 

d c  =  ndeP»  (4.4) 


then  k  =  /2  and  J  is  the  creep  function  in  pure  shear. 

For  our  purposes  it  is  more  convenient  to  write  equation  (4.1)  in  the  expli¬ 
cit  form 


eP  =  /  J(z(t)  -  z(t'))  Sfr  df 
0 

for  reasons  that  will  become  apparent. 


(4.5) 


4- . 1 .  Monotonic  Creep  in  the  Presence  of  a  Constant  Stress  History 


In  this  specific  case 

s ( t )  =  sQ  H { t)  (4.6) 

where  H(t)  is  the  unit  step  function  whose  "derivative"  is  the  Dirac  delta  func¬ 
tion.  In  this  instance  substitution  of  equation  (4.6)  in  equation  (4.5)  gives 
the  creep  response  in  the  simple  form 


(4.7) 


where  sQ  is  the  amplitude  of  the  step  function  of  applied  stress.  It  is 
apparent  from  equation  (4.7)  that  if  the  form  of  J(z)  is  known  then  knowledge  of 
2 ( t )  determines  creep  strain  in  terms  of  the  stress  amplitude  sQ.  We  caution 
that  eP  is  not  necessarily  linear  in  s  since  z ( t )  depends  on  s^  as  we  shall 
demonstrate. 

To  this  end  differentiate  equation  (4.7)  with  respect  to  t  and  use  equations 
(4.2)  and  (4.3)  to  find  that  under  monotonic  conditions 

g(?)  =  ksQ  J'  (z)f'1  (4.8) 

where  J'(z)  is  the  derivative  of  J  with  respect  to  z.  Thus 

C  =  g_1  {k  sQ  J'(z)  f"1  }  (4.9) 

But  from  equations  (4.7)  and  (4.3) 


4  =  k  sQ  J'(z)  z 

Thus  from  equations  (4.9)  and  (4.10) 


(4.10) 


g  1{ksQJ' (z)f"1} 

z  iTyrm 

Equation  (4.11)  gives  z(t)  by  numerical  integration  if  f  is  known. 


(4.11) 


To  assign  analytical  forms  to  the  functions  J(z)  and  g(?)  we  appeal  to 


experiment  and  the  underlying  assumptions  of  endochronic  plasticity.  It  is  com¬ 
mon  experience  that  metals  become  more  strain  rate  sensitive  as  the  temperature 


rises.  However,  the  spirit  of  the  endochronic  theory  is  that  this  change  is 
brought  about  not  by  a  change  in  the  form  of  J(z)  but  by  virtue  of  g(c)  which 
is  evidently  dependent  on  temperature  even  though  this  dependence  is  suppressed 
in  equation  (4.2). 


At  room  temperature  where  rate  effects  are  not  significant  J ( z )  is  repre¬ 


sented  very  closely  by  the  analytical  impression 

J(z)  =  JQ  z°  (4.12) 

where  a  is  the  vicinity  of  0.85  for  a  number  of  metals.  This  form  is  retained 
by  virtue  of  the  above  argument,  at  higher  temperatures. 

Experiments  also  indicate  that  under  monotonic  creep  conditions 

eP  =  F(sq)  t6  (4.13) 

i.e.,  that  the  stress  and  time  dependence  of  creep  strain  are  factorable  and 
that  the  time  dependence  is  represented  very  closely  by  a  power  law.  If  during 
monotone  creep  f  is  a  constant  -  which  was  found  to  be  so  for  one  component  of 
the  creep  -  then  for  equation  (4.13)  to  hold  g(£)  must  also  be  a  power  function 
Thus  we  have  set 


•m 


g(5)  =  c" 

In  view  of  these  stipulations  equation  (4.11)  now  becomes 


(4.14) 


1 


1-m  (a-l)(l-m) 


.m  .  -i  \  m 

f  z  =  (ka  SqJq)  z 


m 


(4.15) 


Special  solutions  to  equation  (4,15).  We  proceed  to  give  some  special  solution 
to  equation  (4.15)  when  (a)  f  is  constant  and  (b)  when  f  is  a  power  function  of 
z.  When  f  is  constant  the  solution  is  given  by  equation  (4.16): 

f™  z  =  A  (k  a  sQ  Jq)  n  t  6  (4.16) 

where 

A  =  (|-r)8  »  S'  =  m/l  +  am  -  a ,  n'  =  6'  (4.17) 

p  m 

In  the  case  where  f  =  fQ  z^,  where  $  >  0,  the  solution  to  equation  (4.15)  is 
given  by  equation  (4.16)  as  before  except  that  now  the  constant  B'  is  given  by 
equation  (4.18): 


(4.18) 


V- 


8 '  =  tn/l  +  am  -  a  + 


and  the  constant  f  on  the  left  hand  side  of  equation  (4.16)  is  now  denoted  by  fQ. 


Knowing  z(t)  one  may  now  calculate  the  creep  strain  e^  by  use  of  equations 


(4.7)  and  (4.12),  in  conjunction  with  equation  (4.16).  Thus 

8 


,ct  m  .n-1  ,  ,  <n  8 

A  fQ  (ka)  ( t 


(4.19) 


where 


1  +  an1  =  n,  8 'a 


8 


(4.20a  ,b) 


Pi scussion.  So  far  we  have  represented  the  creep  strain  by  a  single  integral. 
We  have  also  represented  the  creep  function  J(z),  the  strain  rate  sensitivity 
function  g(?)  and  the  hardening  function  f (?)  by  analytical  forms  of  the  power 
type.  By  analysis  we  then  arrived  at  equation  (4.19)  which  is  basically  of  the 
form 

'P  =  c  (4.21) 


B  sQ"  t6 


where  B  is  a  constant,  whereby  the  monotonic  creep  strain  depends  multiplica- 
tively  on  the  stress  amplitude  to  the  power  n  and  the  time  to  the  power  8.  This 
form  has  appeared  frequently  in  the  literature  where  it  has  been  arrived  at  by 
analysis  of  the  data.  It  does  not  for  all  creep  data  and  certainly  not  over  the 
entire  range  of  stress. 

What  is  important,  however,  is  that  the  creep  strain  depends  on  time 
according  to  a  power  law  (equation  (4.21)  in  accordance  with  observation  as  per 
equation  (4.13)  while  the  dependence  on  stress  is  of  a  more  general  type.  One 
can  change  the  dependence  of  eP  on  Sq  by  changing  the  analytical  form  of  J(z)  or 
g(S)  or  both  but  it  seems  that  this  would  vitiate  the  dependence  of  eP  on  a 
power  function  of  t.  Two  other  avenues  are,  however,  available.  One  is  to 
introduce  a  spectrum  of  intrinsic  times,  as  discussed  previously,  i.e.,  a  series 


of  integrals  on  the  right  hand  side  of  equation  (4.1).  The  other  is  to  intro¬ 
duce  a  stress  dependence  in  the  hardening  function  f  .  This  has  been  found  to 
be  the  case  in  other  materials  such  as  polymers. 

Specifically,  if  one  sets 

f  =  Wz*  (4.22) 

then  under  monotonic  creep  conditions  in  the  presence  of  constant  stress 
equation  (4.19)  will  have  the  form  given  by  equation  (4.13),  for  an  appropriate 
choice  of  the  function  f. 

This  approach  alone,  however,  has  been  found  inadequate  to  describe  creep 
under  cyclic  piece-wise  constant  stress  histories.  This  question  as  well  as  a 
constitutive  equation  involving  more  than  one  intrinsic  time  will  be  discussed 
in  the  next  section. 
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5.  Specific  Consti tui tive  Equations  for  304  Stainless  Steel 


In  the  application  of  the  theory  to  304  stainless  steel  at  600°C  and  speci¬ 


fically  to  the  experimental  data  generated  by  Ohno  et  als.  (Ref.  [4]),  two 


terms  were  retained  on  the  right  hand  side  of  equation  (2.28),  i.e., 


eP  =  J] (z]  )*ds  +  J2(z2)*ds 


(5.1) 


For  the  purposes  of  analysis  and  presentation  of  the  results  it  is  more  con¬ 


venient  to  write  equation  (5.1)  in  the  form 


P  P  P 
e P=  e?  +  e2 


(5.2) 


where 


ep  =  J  (z  )*ds 
r  r  r 


(5.3) 


r  =  1,  2.  In  this  case  two  hardening  functions  exist  in  the  sense  of  equation 


(5.4) 


.  dz  1 

dz  =  7-  *  — 

r  fr  4m 


(5.4) 


where  m  is  a  material  constant  found  to  be  equal  to  0.12.  Also  two  creep  func¬ 


tions  J ^ ( z )  and  J 2 ( z )  are  needed  and  these  were  given  the  analytical  forms  shown 


in  equation  (5.5): 


„  an  a., 

Jl  =  J®  z  ,  J2  =  J®  z 


(5.5  a ,b) 


where  a  =  0.836,  a2  =  1.  =  2.34  x  10  ^  MPa,  J2  =  1.58  x  10  ^  MPa.  It  still 


remains  to  determine  the  form  of  the  hardening  functions.  In  order  to  match  the 


monotonic  data  f^  was  represented  by  a  power  function  of  the  form 


f  2  =  2 


(5.6) 


where  $2  =  0.196. 


On  the  other  hand  could  not  be  represented  as  a  state  function  of,  say,  s 
and  z,  or  any  other  variable  for  that  matter.  Rather,  the  experimental  cyclic 
results  of  Ref.  [4],  gave  strong  indication  that  f^  should  be  given  in  differen¬ 
tial  form  of  the  type 

d  log  f1  =  dF ( | s  | ,  z] ) |z  (5.7) 

1 

Note  that  the  right  hand  side  of  equation  (5.7)  is  not  an  exact  differential  and 
hence  f  is  a  function  of  the  stress  history.  The  physical  implication  of 
equation  (5.7)  is  that  a  change  in  z-j  does  not  affect  f^  if  during  the  change 
the  absolute  value  of  the  stress  s  remains  constant.  A  mathematically  more 
explicit  form  for  f^  is 

d  log  fi  =  (  M,  Zj)  d|s|  (5.8) 

The  logarithmic  form  is  not  fortuitous  but  is  a  consequence  of  the  physics  of 
deformation  kinetics  and  specifically  equation  (3.25)  in  view  of  which 

log  f]  =  M  cj  (5.9) 

d  log  f1  =  8  d  (A  e^)  (5.10) 


The  implication  is  that  in  mechanism  1  the  mean  barrier  height  will  change  when 
the  absolute  value  of  stress  changes  but  not  otherwise.  The  consti tui tive 
description  of  the  material  is  complete  once  the  function  F(|s|,  z^ )  is  known. 
The  function  F  is  given  below  for  various  values  of  |s|  (in  MPa): 


F ( 1 34 . 2 ,  z})  =  2.25  +  1.3(1  -  e'30zl ) 
F  ( 1 20,  z-j )  =  F(  1 34. 2,  z] ) 

F ( 90 ,  z 1 )  =  2.19  +  0.9(1  -  e‘50zl) 
F(60,  z})  =  1.86  +  44(1  -  e"85Z1) 


(5. 1  la ,b,c,d) 


v  .’iV 


V.*. 


2 


In  addition  the  value  of  the  hardening  function  at  zero  stress  and  zero  value  of 
z  is  set  at  0.1054.  This  value  together  with  the  relation. 


log  f ( 0 ,  | s | ^ )  -  log  f ( 0 ,  | s ( 2 )  =  F ( 0 , 1 s 1 1  -  F ( 0 , | s 1 2 ) 


(5.12) 


where  s^  and  s^  are  any  two  stress  levels  determines  f  for  various  values  of  the 
initial  stress  applied  at  the  onset  of  a  creep  experiment. 


*»»  '»»  ni. 


srV 
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6.  Comparison  with  Piece-wise  Constant  Stress  Experimetns  of  Ohno  et  als. 


Monotonic  Creep  Experiments 

The  strain  response  to  a  constant  stress  history  is  obtained  by  application 
of  equation  (4.19)  in  conjunction  with  equation  (5.2).  Specifically 


=  l  A 

r=l 


A 


Or 


(ka  ) 
r 


nr~l 


(Vor’ 


(6.1) 


in  the  presence  of  the  constraint 


6 1  =  S2  =  6*  ni  =  n2  =  n*  mi 


m^  =  m 


so  that  equation  (6.1)  becomes 


P  ,  ,n  .0  r  ,  r  m  n-1  ,  n 
v  =  (s«)  t  I  A  (fAJ  ( ka^)  J, 


r=l 


Or 


Or 


(6.2) 


(6.3) 


In  the  case  of  the  linear  model  (2),  fQ2  is  a  constant.  However,  in  the 
non-linear  model  (1),  fQ1  is  a  function  of  |s|.  This  dependence  is  determined 
by  adjusting  fg,  for  various  values  of  |s|,  so  as  to  obtain  optimal  agreement 
between  theory  and  experiment  in  the  case  of  monotonic  creep.  The  function 
fQ1  (  | s  | )  is  shown  in  Figure  3.  With  all  the  other  constants  known,  the 
descriptive  capability  of  the  theory  is  shown  in  Figure  4. 

Cyclic  Creep  Experiments 

Following  Section  2  let  eP  be  the  r'th  partial  shear  strain  such  that 


e  =  J  (z  )*ds 
r  r  r 


and 


_  m 
eP  =  I  e 
r=l 


(6.4) 


(6.5) 


where  in  our  case  m  =  2.  For  our  purposes  it  is  more  convenient  to  write  equation 
(6.4)  in  the  form 


(a.* ; 


.■'.1*  a.  la-1  a 
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ej  =  /  -  z(tr)]  dt 


r  r 


(6.6) 


in  the  specific  case  of  piece-wise  constant  stress  histories  of  the  type  con¬ 
sidered  by  Ohno  et  als. 


||  =  sQ{«(t)  -  a  Mt-t^  +  a  6 (t-t2)  ....  } 


(6.7) 


where  <S(t)  is  the  Dirac  delta-function,  s^  is  the  initial  stress  amplitude  and  a 


is  a  constant.  In  this  set  of  experiments  two  parameters  sQ  and  "a"  define  the 


history  of  stress  -  in  addition  to  the  reversal  times  t^,  t^  ...  t  . 
Substitution  of  equation  (6.7)  in  equation  (6.6)  and  integration  gives  the 
explicit  result 


«?  -  s0  ivv  +  a  1  M,r  J-(2--o) 


n=l 


r  r  rn 


(6.8) 


or 


*  s0  y(zr> 


(6.9) 


where  y(z^)  represents  the  bracket  on  the  right  hand  side  of  equation  (6.8).  We 


differentiate  equation  (6.9)  to  obtain 


IP 


e'  =  s„  y ( z  )  z 
r  0  r  r 


(6.10) 


_  dy 


where  y  =  .  Use  of  equation  (5.4)  then  gives  the  result 


•p  c 

=  $0  y  f,zm 


(6.11) 


Now  we  take  absolute  values  of  both  sides  of  equation  (6.11)  and  use  equation 
(4.3)  to  obtain 
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*r  «  =  sQ  k  |y  I 

But  in  view  of  equation  (5.4) 

1 


C  =  (f  Z  ) 

r  r 


(6.12) 


(6.13) 


Equations  (6.12)  and  (6.13)  combine  to  give  the  following  differential  relation 

between  dz  and  dt 

r  1 
f  m 


r  dzr 
h  *  .I  -m 

'kso  yhr 


(6.14) 


Integration  of  equation  (6.14)  gives  the  relation  between  zr  and  t. 

Substitution  of  z  (t)  in  equation  (6.8)  then  gives  the  desired  relation  ep(t) 
r  r 

and,  therefore,  eP(t)  upon  use  of  equation  (6.5). 

In  our  particular  case 


Jr  =  °0r  zr 


(6.15) 


Of]  N  ct 

yr  =  “r  JOr  Zr  +  ^  (_1)  Jflr  (3r  '  9rn> 

r  r  Ur  r  i  Ur  r  rn 


(6.16) 


In  Figures  5,  6  and  7  we  show  the  experimental  values  of  eP(t)  obtained  in  Ref. 
[4]  for  the  stress  histories  shown.  Also  shown  are  the  analytical  functions 
eP(t)  obtained  (a)  by  the  use  of  the  present  theory  and  (b)  as  reported  by  Ohno 
et  als.  in  Ref.  [4]  using  their  own  theory. 

The  most  significant  difference  in  the  predictive  capability  of  the  two 
theories  lies  in  their  depiction  of  the  creep  recovery  slope  at  points  of  stress 
reversal.  The  endochronic  theory  predicts  an  infinite  slope,  in  agreement  with 
experiment,  while  the  theory  by  Ohno  et  als.  depicts  a  finite  much  shallower 
slope.  Overall  the  predictive  capability  of  the  endochronic  theory  is  very  good 
for  the  type  of  stress  histories  discussed  here. 


References 


1.  Valanis,  K.C.,  "Endochronic  Theory  with  Proper  Hysteresis  Closure 
Properties,"  Systems,  Science  and  Software,  LaJolla,  CA,  EPRI  Report 
NP-1388,  1980. 

2.  Valanis,  K.C.,  and  Lee,  C.F.,  "Some  Recent  Developments  of  the 

Endochronic  Theory  with  Applications,"  Nuclear  Engineering  and  Design, 
69,  pp.  327-344,  (1982). 

3.  Valanis,  K.C.,  and  Lee,  C.F.,  "Endochronic  Theory  of  Cyclic  Plasticity 
with  -'lications,"  J.  of  Appl.  Mech.  Trans.  ASME  (to  appear). 

4.  Ohno,  J.,  Murakami,  S.  and  Veno,  T.,  "A  Constitutive  Model  of  Creep 
Describing  Creep  Recovery  and  Material  Softening  Caused  by  Stress 
Reversals,"  Journal  of  Eng.  Mat.  and  Technol.,  107,  1,  (1985). 

5.  Chaboche,  J.L.,  Dang  Van,  K.,  and  Cordier,  K.,  "Model ization  of  the 

Strain  Memory  Effect  on  the  Cyclic  Hardening  of  316  Stainless  Steel," 
Trans.  5th  SMIRT,  Lll/3,  Aug.  1979,  Berlin,  Germany. 

6.  Krieg,  R.D.,  Swearengen,  J.C.,  and  Rohde,  R.W.,  "A  Physically-Based 
Internal  Variable  Model  for  Rate— Dependent  Plasticity,"  in  Inelastic 
Behavior  of  Pressure  Vessels  and  Piping  Components,  ed.  T.Y.  Chang  and 
E.  Krempl ,  ASME  PVP-028,  pp.  15-28,  (1978). 

7.  Malvern,  L.E.,  "The  Propagation  of  Longitudinal  Waves  of  Plastic 

Deformation  in  a  Bar  of  Material  Exhibiting  a  Strain-Rate  Effect,"  J. 
Appl.  Mech.,  18,  pp.  203-208,  (1951). 

8.  Haisler,  W.,  "Application  of  an  Uncoupled  Elastic — Plastic  Creek 

Constitutive  Model  to  Metals  at  High  Temperatures,"  in  NASA  Symposium 
on  Nonlinear  Constitutive  Relations  for  High  Temperature  Applications, 
May  19-20,  (1982).  The  University  of  Akron,  Akron,  Ohio,  pp.  185-189. 

9.  Bradley,  W.L.,  and  Yuen,  S.,  "A  New  Coupled  Viscoplastic  Constitutive 
Model,"  in  NASA  Symposium  on  Nonlinear  Constitutive  Relations  for  High 
Temperature  Applications,  May  19-20,  (1982).  The  University  of  Akron, 
Akron,  Ohio,  pp.  217-233. 

10.  Cescotto,  S.,  and  Leckie,  F.A.,  "Determination  of  Unified  Constitutive 
Equations  for  Metals  at  High  Temperature,"  Proceedings  Int.  Conf.  on 
Constitutive  Laws  for  Engineering  Materials,  ed.  C.S.  Desai  and  R.H. 
Gallagher,  Jan.  10-14,  (1983),  Tucson,  Arizona,  pp.  105-111. 


11.  Krempl,  E.,  "The  Role  of  Servocontrolled  Testing  in  the  Development  of 
the  Theory  of  Viscoplasticity  Based  on  Total  Strain  and  Overstress," 
ASTM  STP  765,  pp.  5-28,  (1982). 

12.  Walker,  K.P.,  "Research  and  Development  Program  for  Nonlinear 

Structural  Modeling  with  Advanced  Time — Temperature  Dependent 

Constitutive  Relationships,"  NASA  CR-165533,  Nov.  1981. 

13.  Miller,  A.K.,  "An  Inelastic  Constitutive  Model  for  Monotonic  Cyclic, 
and  Creep  Deformation:  Part  I— Equations  Development  and  Analytical 
Procedures,  Part  II — Application  to  Type  304  Stainless  Steel,"  J.  of 
Eng.  Mat.  Techn.  Trans.  ASME,  98^,  pp.  97-133,  (1976). 

14.  Hart,  E.W.,  "Constitutive  Relations  for  the  Nonelastic  Deformation  of 
Metals,"  J.  Engn.  Mat.  Techn.  Trans.  ASME,  98,  pp.  193-202,  (1976). 

15.  Ohashi,  V.,  Ohno,  N.  and  Kawai,  K.,  "Evaluation  of  Creep  Constitutive 
Equations  for  Type  304  Stainless  Steel  Under  Repeated  Multiaxial 
Loading,"  J.  Eng.  Mat.  and  Technol.,  104,  159,  (1982). 

16.  Murakami,  S.  and  Ohno,  N.,  "A  Constitutive  Equation  of  Creep  Based  on 
the  Concept  of  a  Creep  Hardening  Surface,"  Int.  J.  Solids  and  Struct., 
J8,  597,  (1982). 

17.  Krausz,  A.S.  and  Faucher,  B.,  "A  Kinetics  Approach  to  the  Derivation 
and  Measurement  of  the  Constitutive  Equations  of  Time  Dependent 
Deformation,"  STP  765,  ASTM  (1982). 

18.  Wu,  H.C.  and  Yip,  M.C.,  "Strain  Rate  and  Strain  Rate  History  Effects 
on  the  Dynamic  Behavior  of  Metallic  Materials,"  Int.  J.  Solids  and 
Structures,  26»  PP-  515-536,  (1980). 

19.  Wu,  H.C.  and  Yip,  M.C.,  "Endochronic  Description  of  Cyclic  Hardening 
Behavior  for  Metallic  Materials,"  J.  Appl.  Mech.  Trans.  ASME,  pp. 
212-217,  (1981). 

20.  Lin,  H.C.  and  Wu,  H.C.,  "On  the  Rate-Dependent  Endochronic  Theory  of 
Viscoplasticity  and  its  Application  to  Plastic-Wave  Propagation,"  Int. 
J.  Solids  and  Structures,  T9,  PP-  587-599,  (1983). 

22.  Valanis,  K.C.,  "Proper  Tensorial  Formulation  of  the  Internal  Variable 
Theory,"  Archives  of  Mechanics,  29,  173,  (1977). 

23.  Valanis,  K.C.,  "Fundamental  Consequences  of  a  New  Intrinsic  Time 
Measure:  Plasticity  as  a  Limit  of  the  Endochronic  Theory,"  Arch,  of 
Mech.,  32,  pp.  171-191,  (1980). 

24.  Kraus,  A.S.  and  Eyring,  H.,  "Deformation  Kinetics,"  Wiley,  N.Y. 
(1975). 

25.  Valanis,  K.C.  and  Lalwani,  S.,  "Thermodynamics  of  Internal  Variables 
in  the  Content  of  the  Absolute  Reaction  Rate  Theory,"  J.  Chem.  Phys., 
67,  3980,  (1977). 


Fig.  2  Energy  Distribution  in  Vicinity  of  e0 

€;  =  ek*  ( 

wr-  ( v  US 

Wf  -  (i.-  if)  8 


Pf‘Pb  = 


2  e 


/*A« 


77W 


1  -e 


sinh  /3fP 


Fig.  4  Monotonic  Creep  Tests 


Strain  Response  to  Piece-wise  Constant 
Stress  History  (/5  S  =  120  Mpa,  a  =  1.0) 


Strain  Response  to  Piece-wise  Constant 
Stress  History  (/ 5  S  —  1 20  Mpa,  a  =  1.5) 


Present  Model 


Response  to  Piece-wise  Constant 


MAILING  LIST 


NAME 

STREET  ADDRESS  AND  AFFILIATION 

CITY 

STATE 

ZIPC0DE 

Dr.  Charles  D.  Babcock 

2392  Mors! ay  Road 

Altadena 

CA 

91001 

Dr.  Paul  Cernocky 

Shell  Development  Co.,  PO  Box  481 

Houston 

TX 

30332 

Dr.  Peter  C  T  Chen 

Watervliet  Arsenal,  US  Army  R&D 
Comm  LCW 

Watervl iet 

NY 

12189 

Dr.  R 

M  Christensen 

Lawrence  Livermore  Lab,  PO  Box  808 
(L-338) 

Livermore 

CA 

94550 

Prof. 

Y  F  Dafalias 

Dept,  of  Civil  Engng,  Univ.  of 

Cal i fornia-Davi s 

Davis 

CA 

95616 

Dr.  L 

W  Davison 

Sol  id  Mechs  Dept  1 530, 

Sandia  Natl  Labs 

A1 buquerque 

NM 

87185 

Dr.  John  K.  Dienes 

Theoret  Div. ,  MS  216, 

Los  Alamos  Scient  Lab 

Los  Alamos 

NM 

87545 

Dr.  Stanley  Fistedis 

Reactor  Analy  &  Safety  Div., 
Argonne  Nat'l  Lab 

Argonne 

IL 

60439 

Prof. 

J  E  Fitzgerald 

Sch  CE,  Georgia  Tech 

Atlanta 

GA 

30332 

Prof. 

M  E  Gurtin 

Mathematics,  Carnegie-Mellon 

Univ. 

Pittsburgh 

PA 

15213 

Prof. 

I  eonard  Herrmann 

CE  Dept.  Bainer  Hall , 

Univ.  of  Cal ifornia 

Davis 

CA 

95616 

Dr.  Walter  Herrmann 

Org  1500,  Sandia  Nat'l  Labs 

A1 buquerque 

NM 

87185 

Prof. 

T  J  R  Hughes 

Applied  Mechanics,  Duran  281, 
Stanford  Univ. 

Stanford 

CA 

94305 

Melvin  F.  Kanninen 

Eng  &  Matls,  Southwest  Res  Inst, 
6220  Culebra  Road 

San  Antonio 

TX 

78284 

Prof. 

Dusan  Krajcinovic 

Civil  Engr.  P  0  Box  4348 

Univ  of  Ill. 

Chicago 

IL 

60680 

Prof. 

F  A  Leckie 

T&AM,  296  Tlabot  Lab 

Univ  of  Ill,  104  S.  Wright 

Urbana 

IL 

61801 

Dr.  J 

K  Lee 

Engr  Mechs,  Ohio  State  Univ 

Columbus 

OH 

43210 

Prof. 

Frederick  F.  Ling 

ME  AE  &  Mech,  RPI 

Troy 

NY 

12181 

Mr.  J 

Lubl iner 

CE  Dept.  Univ  of  California 

Berkel ey 

CA 

94720 

Prof. 

L  E  Malvern 

Engr  Sc i  Dept,  231  Aero  Bldg, 

Univ  of  Florida 

Gainesvil 1 e 

FL 

32611 

Prof. 

X  Markenscoff 

Mech  &  Environ  Engng, 

Univ.  of  California 

Santa  Barbara 

CA 

93106 

Prof. 

S  Nemat-Nasser 

Civil  Engng,  Northwestern  Univ. 

Evanston 

IL 

60201 

Prof. 

J  Tinsley  Oden 

Ticom  WRW  305,  Univ  of  Texas 
<?  Austin 

Austin 

TX 

78712 

Dr.  Steven  Peng 

Jet  Propulsion  Lab, 

67/201  Caltech  Bldg. 

Pasadena 

CA 

91103 

NAME 


STREET  ADDRESS  AND  AFFILIATION 


CITY 


STATE 


ZIPCODE 


«• 

Dr.  Jerry  T  Pendera 

Inst  for  Experimental  Mechs, 

Univ.  of  Waterloo 

Waterloo,  Ont 

CANADA 

N2L  3G1 

b 

Prof.  Karl  S.  Pister 

CE  Dept,  Univ  of  California 

Berkeley 

CA 

94720 

J  W  Rudnicki 

Civil  Engng,  Northwestern  Univ 

Evanston 

IL 

60201 

Prof.  Jerome  L  Sackman 

Civil  Engng,  Univ.  of  California 

Berkeley 

CA 

94720 

r 

Dr.  Anthony  San-Miguel 

2775  Mesa  Verde  Dr  East,  Apt  P-212 

Costa  Mesa 

CA 

92626 

P 

Prof.  S  S  Sternstein 

Mata  Engr  Dept.  RPI 

Troy 

NY 

12181 

Mm 

Prof.  Gerald  A  Wempner 

Engng  Sci  &  Mechanics, 

Georgia  Tech 

Atlanta 

GA 

30332 

1 

Prof.  Kaspar  J.  William 

Campus  Box  428,  Univ  of  Colorado 

Boulder 

CO 

80309 

Mr.  Lawrence  K  Yu 

Inti  Paper  Co.  Box  797 

Tuxedo  Park 

NY 

10987 

_ LXCI.AS.Sm--'') _ _ _ . _ _ _ 

SECURITY  CLASSIFICATION  OF  TH'S  PAGE  rHTi.n  r>«tf.  hntered) 


REPORT  DOCUMENTATION  PAGE 


m  172  70  S' 


I.  REPORT  NUMBER 


4.  TITLE  (end  Subtitle) 

An  Endochronic  Rate-Sensitive 
Constitutive  Equation  for  Metals 
Application  to  Generalized  Creep 


7.  Au  T HOR^J 

K.  C.  Valanis,  C.  F.  Lee,  and  S.  D.  Lee 


MASTER  COPY  -  FOR  REPRODUCTION  PURPOSES 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT’S  CAT  ALOG  NUMBER 


s.  type  of  report  a  period  covered 

Technical  Report 
12-1-84  till  11-30-85 


6.  PERFORMING  ORG.  REPORT  NUMBER 


8.  CONTRACT  OR  GRANT  NUMBER(m) 


DAAG  29-85-K001*? 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 


University  of  Cincinnati 
Cincinnati,  Ohio  45221 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  &  WORK  UNIT  NUMBERS 


It.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

U.  S.  Army  Research  Office 
Post  Office  Box  12211 


12.  REPORT  DATE 

August  1986 


13.  NUMBER  OF  PAGES 


14.  MONITORING  AGENCY  NAME  8  ADDRESS^//  different  from  Controlling  Office)  15.  SECURITY  CLASS,  (of  thle  report) 


U.S.  Army  Research  Office 
P.0.  Box  12211 

Research  Triangle  Park,  NC  27709 


16.  DISTRIBUTION  STATEMENT  (of  l/il.  Report) 


Unclassified 


15a.  DECLASSIFICATION/  DOWN  GRADING 
SCHEDULE 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  mbetrect  entered  In  Block  20,  It  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 

The  view,  opinions,  and/or  findings  contained  in  this  report  are 
those  of  the  author(s)  and  should  not  be  construed  as  an  official 
Department  of  the  Army  position,  policy,  or  decision,  unless  so 

_  prmf  P.d  bv  nrhpr  rlnrumpntatinn. 


19.  KEY  WORDS  (Continue  on  reveree  eldm  It  neceeeery  end  Identify  by  block  number) 

Endochronic,  plastic,  viscoplastic,  rate-sensitive  constitutive  equation, 
generalized  creep. 


20.  ABSTRACT  CCaattmie  ere  revere*  etde  ft  neceeeery  mad  Identity  by  block  number) 

A  constitutive  equation  is  proposed  with  a  view  to  describing  the  rate 
dependent  mechanical  response  of  metals  at  high  temperatures.  The  equation 
is  of  the  endochronic  type  and  derives  its  physical  foundations  from 
deformation  kinetics.  Of  importance  is  the  fact  that  hardening  is 
associated  with  a  change  in  the  energy  barriers  brought  about  by  the 
inelastic  deformation  of  a  metal.  The  equation  is  used  to  describe  the 
results,  by  Ohno  and  his  associates,  of  experiments  on  the  creep  response 


EDITION  OF  »  MOV  65  IS  OBSOLETE 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  Dele  Entered) 


SECURITY  CLASSIFICATION  OF  THIS  PAGE(T*7i«n  Data  Entarad )  _ ___ 


of  metals  to  piece-wise  constant  stress  histories.  The  metal  in  the  present 
case  is  304  stainless  steel  at  600°  C.  It  is  shown  that  the  theory  gives 
analytical  results  that  are  in  close  agreement  with  experiment. 

Of  consequence  is  the  fact  that  the  constitutive  equation  applies  to 
three-dimensional  stress  or  strain  histories  and  is  thus  not  limited  to 
those  stress  histories  associated  with  creep. 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  P AGEfWTi»n  Data  Entarad) 


